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Abstract 

The propagation of an external transverse magnetic signal acting 
locally on a Id chain of spins generates a disturbance which runs 
through the system. This quantum effect can be interpreted as a 
classical traveling wave which contains a superposition of a large set 
of frequencies depending on the size of the chain. Its local amplitude 
fixes the size of the z-component of the spins at any location in the 
chain. The average and maximum value of the group velocity are 
determined and compared with the transmission velocity fixed by the 
Lieb- Robinson upper bound inequality. 
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1 Introduction 

Modern communication needs the development of technologies which allow 
fast, faithful and robust information manipulation by means of adapted trans- 
fer devices. These needs generate intense research in different fields of con- 
densed matter physics in order to develop efficient material tools on which 
information can be safely propagated. 

Indeed, the last decade saw the development of new several techniques and 
methods in the field of quantum information technology such as ultra-fast and 
flexible optical manipulations of spins in magnetically ordered materials [H 
121 E], transfer of quantum states and qubits along linear systems and different 
techniques lilElElITllIlEllIOllIIliniinilll], superconducting devices [151 US] , 
transfer of entangled states [T71 [HI [191 ISHl [21], investigations concerning 
echoes and decay times [22], transition times [231 [2l], propagation velocities 
in solid state systems [2Sl[2Zl[2Sl[2a[3nilMl[321[33]. The subject raised 
also more fundamental points concerning the properties related to the non- 
equihbrium behaviour of such systems [Ml [351 [26] . 

The present work uses a simple model, a finite XY chain submitted to 
even boundary conditions in its ground state. It works here as the material 
support for the propagation of external perturbations along the ring chain. 
The perturbations are chosen to be local excitations induced on a single 
spin over a short interval of time and the propagation of the signal in time 
is followed over the whole chain. The interest here lies essentially on the 
estimate of the propagation velocity of a signal through the chain. The per- 
turbation behaves like a traveling wave characterized by a set of frequencies 
which depend on the size and eigenenergies of a strongly correlated system. 
This quantity can be obtained in an analytic form and group velocities are 
defined. Their order of magnitude are obtained by means of numerical ap- 
plications and compared with the upper velocity obtained in the application 
of the Lieb-Robinson (LR) approach to the transmission of information in 
discrete quantum systems. 
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Section 2 introduces the model and the time evolution of the z-component 
of local spins along the chain. In section 3 different group velocities are 
defined and the expression of the LR bound velocity is recalled. Numer- 
ical applications are presented in section 4 and the results concerning the 
propagation are discussed. Section 5 is concerned with some considerations 
concerning the control of the local spin amplitudes during the propagation 
process and the implementation of sequences of local excitations in the sys- 
tem. Conclusions are drawn in section 6. 

2 The physical system 
2.1 The model 

The Hamiltonian of the finite Id quantum chain with sites and asymmetry 
parameter 7 reads 

H, = J/2(l + 7) 5^ atat^^ + J/2(l - 7) E ^'^1 

(i) 

In the sequel even periodic boundary conditions are introduced. The inter- 
action strength is chosen to be J = 1 and h = 1. In Fourier space each site 
[k = 1, N/2] contains four single spin states. The energy spectrum can be 
determined analytically in terms of single particle energies [10| Hi] . Diago- 
nalization of the single particle (4 * 4) Hamiltonian leads to the stationary 
orthonormalized eigenstates 

1^^) = a1\0)k + a',atat,\0)k 

m) = a+jO), (2) 

where \0)k is the spin vacuum, creates a spin at site k, and are 
complex amplitudes which are shown in Appendix A. 
The corresponding eigenvalues read 
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COS (pk - [(cos (pk - hof + 7^ sin^ (pk] 
cos (pk + [(cos (pk - hof + 7^ sin^ cpk] 

J, (4) 

COS (pk = el' 



1/2 



1/2 



(3) 



where (pk = 2nk/N. 

The iVt/i spin of the system experiences an external magnetic field hi 
which induces a perturbation 



with = The external field acts on the z component of the spin 

at site N over a finite time interval fixed by th and starting at t = 0. The 
generated signal is transmitted to the whole chain of connected interacting 
spins. The localized action of an external magnetic field may be an experi- 
mental challenge. There exists however experimental work which shows that 
this might be possible in the not too far future ESI ES]- 

2.2 Time evolution of the individual spins in the chain 

The z component of each spin site n evolves in time as S^{t) 



The physical situation considered here corresponds to an excitation in- 
duced by H[^\t) which relaxes over a time th- This time is chosen as 
small compared to the time unit. In this regime a perturbative treatment of 
the expectation value of S^{t) makes sense. The perturbation is essentially 
governed by the parameter th and hiTn which can be small in practical ap- 
plications if Th is a very short relaxation time even for sizable magnetic fields 
h. 

Using a perturbative treatment up to second order in the perturbation 
development provides an analytic expression of the time evolution of the 
spins along the chain. The ground state |\l/o) of the system is the product of 
the single particle states. A second order perturbation expansion gives the 



i7r^(t) = /iiexp(-t/r^)5^ 



(4) 



^^(t)=exp(i [ H{t')dt')S'^{0)exp{-i [ H{t')dt') 



(5) 




where H{t) = Ho + Hl''>{t). 
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expectation value {S^{t)) of the z-component of the spin with respect to the 
ground state dehvers three contributions 

(sm) - {s:Y'^ + {smY'^ + (smY'^ (6) 

where 

N/2 

(5^)(0)^2/7V^|a,T-l/2 (7) 
fe=i 

N/2 

= 1/2 5] 4'''^[^«Nn«.,) + sin(^-,,)] 

k,l=l 

+Bki[cos{ij+,^) + cos(^„^,,)]] + C{t, th) (8) 

where f 2'^''-' arc transition matrix elements between different single spin states 
{k,l), Aki and Bki amplitudes, C{t,TH) a time-dependent and state indepen- 
dent contribution which relaxes to a constant over a time interval t > th- 
The phases read 

V'n;w = ^^(ki)t + nK{k, I) 

^n;kl = ^(W)^ - ^^(^' (9) 

with U(ki-) = e^^'* -|- e[^^ — e^^^ — e^"* and K{k, I) = {(pk — 4>i)- The expressions 
of the matrix elements ^2'^''^ and the amplitudes Aj^i and B^i are shown in 
Appendix B. 

The expression of the second order contribution {S^^it))^"^^ is given by a 
large number of terms. Present applications show that their contributions 
are sizably smaller than those of the first order when t ^ th- In Appendix 
C we estimate their order of magnitude. 

3 Group velocities 

The time-dependent contribution of {S^{t))^^^ given by equation (8) can be 
interpreted as a sum of classical traveling waves characterized by frequencies 
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co(^ki) and phases (0^ — (pi). Considering the analytic expressions of the co(ki) 
and the energies e^*^ given above one sees that these quantities are incom- 
mensurate. Hence the oscillations will not be periodic in time. The group 
velocity for a fixed wave number K{k, I) characterizes the propagation of the 
waves. This quantity is defined as 

Vgr = duj{K)/dK (10) 

For a fixed value of K{k, I) more than one value of u{k, I) is possible since 
K{k,l) depends on the difference k — I. In the sequel we consider u"^"'^{K) 
which correponds to maximum values of us for a given set of couples {k, I) 
corresponding to a fixed K and the averages < u > over all possible couples 
{k, I) generating the same K. 

The group velocity is a measure of the speed at which the external per- 
turbation propagates through the medium, here in the case of a spin chain. 
By means of Eq.(lO) it is possible to define and obtain an order of magnitude 
of this quantity and an upper value of the propagation of a signal in this type 
of quantum devices. Explicit calculations are presented and discussed below. 

On the other hand the Lieb-Robinson procedure [26] allows the determi- 
nation of an upper limit of the velocity v with which a signal is transmitted 
in a discrete system from an initial point at a fixed initial time to a final 
point at a fixed later time This velocity depends on metric 

factors which define distances between sites of the system and the norm 
of the Hamiltonian matrix. The limit velocity is obtained by means of a 
variational procedure depending on a real positive parameter a. Working 
out the expression and determining the upper bound by varying a leads to 
the expression 

VLR = eN\\H\\/2 (11) 

where e is the exponential constant, N the number of sites of the system. 
Some details of the derivation are given in Appendix D. In section 4.3 we 
work out vlr and compare it with Vgr- 

4 Explicit calculations 
4.1 Preliminaries 

The time evolution of a fixed spin and the behaviour of the spins in the chain 
are worked out for a spin chain of finite length. The physical conditions under 
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Figure 1: z-components of the spins n for t=l. 



which the perturbative treatment discussed above is vahd are reahzed in the 
present apphcations. The zeroth order contribution to the z-componcnt of 
the spin (5'^)''°'* is time- independent and the second order term (S^)^"^^ is 
smaller than (S^)^-^^ by two orders of magnitude. Due to the symmetries 
induced by the boundary conditions imposed on the chain the oscillation 
frequencies u!{k, I) are symmetric with respect to k and I belonging to the 
intervals [1, N/2\, [—1, —N/2]. Hence the calculations can be restricted to the 
interval [1, N/2\. In the sequel we consider a finite chain of length = 100 
and a set of parameters (/iq = 0.5, hi = 1.0,7 = 0-5 which correspond to a 
system which lies above the critical point ho = he = —0.5 after the action of 
the short excitation corresponding to th — 10~^. 

4.2 Spin oscillations along the chain 

The time evolution of 



for fixed time and a set of sites n = 1 to 50 is shown in Figs. 1 and 2. 

For short times the essential part of the excitation resides in the sites with 
low n which are close to = 100 due to the periodic boundary conditions. 
For large times the longitudinal part of all the spins is affected and shows 
a somewhat irregular and non repetitive pattern due to the incommensura- 
bihty of the values of u and K(k, I) discussed above. The calculations have 
been repeated for an excitation hi > above ho — he — —1/2 where the 
system experiences a phase transition as well as an excitation which leads 



(12) 
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Figure 2: z-components of the spins n for t=200. 



to a crossing of the critical point when t goes to infinity. In both cases no 
specific behaviour of the spin component is observed. Individual spins which 
are affected by local excitations are qualitatively insensitive to the global 
effects induced phase transitions which affect the whole system. 

Figs. 3 and 4 show the behaviour of {S^^-^it))'^^'^^ and {S^^^^^it))'^^'^^ over 
an interval of time t = [1,200]. It is seen that the evolution of the spins 
in time corresponds globally to a wave which spreads over all spins and 
decreases necessarily in amplitude since the total magnetization of the system 
is constant and the perturbation is small. 



4.3 Group velocities 

Figs. (5) and (6) show the values of the maximal and average group velocities 
v^^^'-'iK) and < v(^) > S^^d K = [l,N/2 - 1], K{k,l) = \k - l\. The 
parameters of the model arc those given in section 4.1. One observes sizable 
variations of these quantities over the range of K values which shows a large 
dispersive behaviour of the system due to the fact that the frequencies show a 
strongly non- linear dependence on K{k, I). This quantity is is itself a mixture 
of many single particle energies. The absolute value of the maxima are 1.5 for 
Vg^"'^{K) and 0.55 for < Vgr{K) > which correspond to a maximum transit 
time between two sites of 0.66 time units for Vg^""^{K) and 1.8 for < Vgr{K) >. 
These quantities are independent of the amplitude of hi, the velocities depend 
on the characteristics of the system only. 

Using the expression of the Lieb-Robinson velocity for a signal which 
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evolves over the whole system one finds vlr — 1.02 * 10"^ (iV — 100 and 
\\H\\ = 75) in the same physical units as those used for Vyj.. It appears that 
there is a large quantitative difference between the two velocities. One can 
make the following comments. 

The group velocity Vgr generated by an external disturbance is a quantity 
which is local in K space and depends on single particle energy differences, 
see section 2.2. On the other hand v^^ is the maximum speed with which a 
signal emitted at some time at some place is transmitted through the system 
and observed at some other place at a later time. It depends explicitly on the 
extensive properties of the system, the maximum of its total energy (||-f^||) 
and its size (A^). The transmission between two neighbouring sites {N — 2) 
one would lead to tij^/j = 4.07. 

The expression of Vlr depends on the geometric stucture of the system 
only through its total energy which although different should be of the 

same order of magnitude for open and closed systems. This effect is certainly 
weaker than the difference between Vgr in a closed and an open system due 
to the wave nature of the propagation. 

5 Information transmission along the chain 
5.1 Spin amplitude transportation 

Consider the transmission of the amplitude of the z-component {S^{t))^^^ 
located at site N at time t to some other site in the chain. This may be 
realized in the following way. 

At a given initial time t and fixed hi and th one applies an excitation at 
site N which generates a perturbation {Sf^{t)y^\ At time t + At, {Sff_i{t + 
At, /ii))^^^ has a well defined value. If one applies a further excitation h[ at 
site N — 1 and time t + At such that 

(5^_,(t + At,/^;))« (13) 

the amplitude of the spin component at site — 1 at time t + At will be the 
same as the amplitude at site N at time t. The field h[ needed is easily fixed 
since any z-component of a spin {Sn {t + At, hi)) depends linearly on the 



9 




-0.1 - 



50 100 150 200 

t 

Figure 3: Time dependence of the z-component {S^{t))^^'^^ of the spin at site 
n=l. 

field hi under the physical conditions described above. The operation can be 
iterated over a finite interval of sites of a given length. 

The fidelity with which these operations can be realized as presented here 
are of course subject to limitations. Indeed the excitation must be such that 
higher than first order corrections are effectively negligible. Furthermore the 
physical realization of a transmission device of this type is also subject to the 
precision with which the operations described above can be realized. This 
concerns essentially the mechanical and thermal fluctuations of the external 
field applied at site N. 

5.2 Finite sequences of excitations 

In practical applications the excitation to be transmitted from site to some 
site m through the system might be repeated a finite number of times starting 
from an initial time t = 0. Consider the case where a set of n excitations are 
sent from site N over regular intervals of time to ^ th- For times t > nto 
there are two contributions to {S^{t)y^\ (S^it))^^^ = {SL{t)Y^''^ + {S^{t)Y^^^ 
which read 
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t 

Figure 4: Time dependence of the z-component of the spin {S^{t))^^'^^ at site 
n=40. 



(1 - cos(a;(w)(t - pto)) exp(-(i - pto))) 

sin(c<;(fci)(t -pto))exp(-(t - pto))] (14) 
where D{k, I) — [l/r^ + ^fki)] and 



2^. 2 

k,i 



with 



Sn{k, l)[Aki[{sin{ipl^.f^i) + sin((^„_^.fc;)] + 

+C«[(cos((^J^.^^) + cos((^-^.^^)]] (15) 



^ ^ sin[a;(fc,)/2to] ^ ^ 



and 
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vtm;ki = ^{ki){t - nto/2) + mK{k, I) 

fn^m-M = ^(ki)it - nto/2) - mK{k, I) (17) 

The observable {S^it))^^"^ is known at each time t and can in principle 
be renormalized to a fixed value as discussed above in the case of a single 
excitation. 

6 Summary and outlook 

The propagation of a local perturbation generated on a quantum spin chain 
at temperature T = has been investigated. The perturbation generated by 
an external magnetic field starts at an initial time t = at some local spin 
site of a Id chain and affects the transverse component of the other spins. 
The conditions imposed on the relaxation time of the external field allow a 
perturbative treatment. 

The action of the external field induces time-dependent oscillations of 
the transverse components along the spin sites of the chain which behaves 
like a ring due to the imposed boundary conditions. The oscillations of the 
spin components show the complex behaviour of classical traveling waves 
with a large number of harmonic modes due to the spectral properties of the 
physical system on which the perturbation propagates. The case where finite 
sequences of perturbations are locally emitted at fixed time intervals of equal 
length has been worked out. 

In order to characterize the propagation of the perturbation different 
group velocities have been introduced. These velocities show a strong dis- 
persive behaviour far from sharp wave packets. They reflect the complex 
nature of the strongly coupled spins in the chain. It comes out that these 
velocities are sizably smaller than the Lieb-Robinson propagation velocity 
bound deflned in [281 129]. Reasons for this have been presented and dis- 
cussed in section 4. 

The chain undergoes a phase transition for a speciflc value he of the 
internal magnetic fleld which acts on the transverse components of the spins. 
The propagating wave is not qualitatively influenced by this property. It 
is due to the fact that the propagation concerns individual spins and the 
external fleld does not affect the global behaviour of the system. This point 
is analysed in section 4. 
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Figure 5: Maximum group velocity as a function of K(k,l) for a chain with 
N=100 and the parameters given in the text. 

The way in which the value of the transverse component of a spin at 
one site can be retrieved at another site has been investigated under ideal 
noiseless conditions. A more realistic approach should include the presence 
of noise which can be induced by initial thermal effects and fluctuations of 
the external magnetic field acting on the system. It would also be of interest 
to consider the case of a chain coupled to a heat bath at temperature T ^ 
at time t = which would introduce the action of excited states of the chain. 

The present study concerns systems which arc influenced by a fleld which 
acts over a small amount of time and cannot produce a reversal of the spins. 
The extension of the present investigation to chains excited by strong flelds 
able to affect sensibly or even to reverse the transverse spin direction would 
request a non-perturbative treatment. Such a case needs further develop- 
ments. 

Finally the choice of even boundary conditions on the end sites of the 
chain closes it in practice. The behaviour of the propagation in an open 
chain may lead to a different behaviour of the perturbation. In the prospect 
of practical applications it would be particularly interesting to investigate 
this point. 

Thanks arc due to J.-Y. Fortin for a careful reading of the manuscript. 

Appendix A: explicit expressions of the sin- 
gle particle amplitudes 
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Figure 6: Average group velocity as a function of K(k,l) for a chain with 
N=100 and the parameters given in the text. 

The energy spectrum of the Id XY chain of length reduces to a product 
of single particle energies with 4 states for a given momentum {k = [1, A^/2]) 
in Fourier space [ID]. The local matrix for fixed k can be analytically diago- 
nalized leading to the eigenvalues given by the expressions shown in Eq.(3). 
The amplitudes of the eigenvectors corresponding to the states 1 and 2 are 
given by 

a1 = (el')-(2cos(0fc-/io))/Ar(i) 

= (-z5,/2)/iV(i) 
/3f = (^5p/2)/iV(2) 

^2 = {ef-h,)/N^,) (18) 
with 5p = —27 sin(0j!c) and the normalization factors 

Ar(,) = [1/45J + (6« - (2COS0, - MT^' 

= + (ef - MT^' (19) 



Appendix B: matrix elements and amplitudes 
entering {sf^\t)) 
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The transition matrix elements V2^'''^ of the first order contributions to 
{Sj^^\t)) are given by 

v^^''^ = {Mk)\v\^^m 

t;f') = (vl>o(fc)|f|vl>,.(A;/)) (20) 

where I'^o^k)), I'^ij^kl)) are the single state eigenf unctions at sites 

(fc,/) in Fourier space in the ground state |0 > and excited states \i > and 
\j >,V the time-independent part of the diagonal and non-diagonal elements 
of the perturbation interaction V{t) = /ii exp(— t/rj7)c^CAr expressed in the 
second quantization formalism in ordinary space [40] . 

The amplitudes Aki and Bki depend on the structure of the chain and the 
applied magnetic field at site N. They read: 



A 



kl 



1/th 



a(m) 



= ^^V^ (21) 

where Al^'-'^je = e^*^ + ep'' — e[.^^ — e[^\ i and j corresponding to excited states. 

Appendix C: second order contributions to 
the tranverse component of spin n 

The order of magnitude of the matrix elements given by the coefficients 
linked to the time-dependent part of the different contributions of second 
order are given by the strength factors Str 

Si = hlrl 
S2 = hi{l/D') 
Ss = hl/{THD) 
D ~ (1/4 + Ae^) (22) 

where Ae ~ e^*) — e'^-'^ (z, j) corresponding to single particle ground or excited 
states. The single particle state energies are of the order of unity. In the 
numerical applications hi is chosen to be of the same order of magnitude and 
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th is four order of magnitude smaller. As a consequence the second order 
contributions to {S'^) are two to three orders of magnitude smaller than the 
first order ones. This justifies their neglect in the numerical estimates. 

Appendix D: estimate of the Lieb- Robinson 
bound for the chain. 

The expression of the velocity reads 

Va = iniCa\\^a\\la (23) 

a 

Here a is a real positive number, Ca depends on the distances between discrete 
sites {x,y,z) of the system and 

= \\H\\/D,{x,y) (24) 

where is the norm of the Hamiltonian matrix and Da{x,y) a measure 
of the distance between sites x and y multiplied by an exponential factor 
exp{-a\x-y\) [271128]. 

Working out this expression leads to the expression 

Va = miN/2exp{aN)^{a = l)/aN (25) 

a 

which shows an extremum for a = 1/N. This extremum is a minimum. The 
final result is given by vlji in section 3. 
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